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Abstract 

Inelastic photoproduction of J/ip particles at high energies is one of the processes 
to determine the gluon distribution in the nucleon. The QCD radiative corrections to 
the color-singlet model of this reaction have recently been calculated. They are large at 
moderate photon energies, but decrease with increasing energies. I compare the cross 
section and the J/ip energy spectrum with the available fixed-target photoproduction 
data. Predictions for the HERA energy range are given which demonstrate the sen- 
sitivity of the result to the parametrization of the gluon distribution in the small-x 
region. 



* Talk presented at the Workshop on "Heavy Quark Physics", Bad Honnef, FRG, Dec. 1994. 
f E-mail: mkraemer@desy.de 



The measurement of the gluon distribution in the nucleon is one of the important goals of 
lepton-nucleon scattering experiments. The classical methods exploit the evolution of the 
nucleon structure functions with the momentum transfer and the size of the longitudinal 
structure function. With rising energies, however, jet physics and the production of heavy 
quark states become important complementary tools. Besides open charm and bottom pro- 
duction, the formation of J/ifj bound states in inelastic photoproduction experiments 

-i + U ^ J/ij + X (1) 

provides an experimentally attractive method since J/ip particles are easy to tag in the 
leptonic decay modes. 

Many channels contribute to the generation of J/ip particles in photoproduction exper- 
iments 10], similar to hadroproduction experiments. However, no satisfactory quantitative 
picture has emerged yet and the production of a large surplus of ip' particles in pp collisions 
awaits the proper understanding. Theoretical interest so far has focussed on two mecha- 
nisms for J /ip photo- and electroproduction, elastic/diffractive [|],|[ and inelastic production 
through photon-gluon-fusion While by the first mechanism one expects to shed light 

on the physical nature of the pomeron, inelastic J/ip production provides information on the 
distribution of gluons in the nucleon 0. The two mechanisms can be separated by measur- 
ing the J/ ip energy spectrum, described by the scaling variable z = p ■ / p ■ with p, k^^^ 
being the momenta of the nucleon and J/ip, 7 particles, respectively. In the nucleon rest 
frame, z is the ratio of the J/ip to the 7 energy, z = E^/E^. For elastic/diffractive events 
z is close to one; a clean sample of inelastic events can be obtained in the range z < 0.9 
0. The production of J/ip particles at large transverse momenta is dominated by gluon 
fragmentation mechanisms 0. Additional production mechanisms, such as BB production 
and the "resolved photon" contributions at HERA can be strongly suppressed by suitable 
cuts 0. 

Inelastic J/ip photoproduction through photon-gluon fusion is described in the color- 
singlet model through the subprocess 

-i + g^J/i^ + g (2) 

shown in Fig.|l]. Color conservation and the Landau- Yang theorem require the emission of a 
gluon in the final state. The cross section is generally calculated in the static approximation 
in which the motion of the charm quarks in the bound state is neglected. In this approxima- 
tion the production amplitude factorizes into the short distance amplitude 7 + (7 — > cc + (yf, 
with cc in the color-singlet state and zero relative velocity of the quarks, and the cc wave 
function v^(0) of the J /if) bound state at the origin which is related to the leptonic width. 
When confronted with photoproduction data of fixed-target experiments , the theoretical 
predictions underestimate the measured cross section in general by more than a factor two, 
depending in detail on the J/ ip energy and the choice of the parameters . The discrepancy 
with cross sections extrapolated from electroproduction data |10;11| is even larger. 



The lowest-order approach to the color-singlet model demands several theoretical refine- 
ments: (i) Higher-order perturbative QCD corrections; (ii) Relativistic corrections due to the 
motion of the charm quarks in the J/ip bound state; and last but not least, (iii) Higher-twist 
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effects which are not strongly suppressed due to the fairly low charm-quark mass. While 
the relativistic corrections have been demonstrated to be under control in the inelastic re- 
gion 



1^, the problem of higher-twist contributions has not been approached so far. The 



calculation of the higher-order perturbative QCD corrections has been performed recently 
13| . Expected a priori and verified subsequently, these corrections dominate the relativistic 



corrections in the inelastic region, being of the order of several as{Mj^^) ~ 0.3. In the first 
step of a systematic expansion, they can therefore be determined in the static approach [|14|. 

A detailed analysis of the 0{aal) corrections to inelastic J/ip photoproduction is the 
subject of a forthcoming publication |T^; first results have been presented in Ref. |]16[. In this 
short note the implications of the higher order QCD corrections for the partonic cross sections 
are discussed and the J/ip energy spectrum is compared with the available fixed-target 
photoproduction data. In addition, predictions for the HERA energy range are given which 
demonstrate the sensitivity of the result to the parametrization of the gluon distribution in 
the small-x region. 

Generic diagrams which build up the cross section in next-to-leading order are depicted 
in Fig.||. Besides the usual self-energy diagrams and vertex corrections for photon and 
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Figure 1: Generic diagrams for inelastic J/ip photoproduction: (a) leading order contribu- 
tion; (b) vertex corrections; (c) box diagrams; (d) splitting of the final state gluon into gluon 
or light quark-antiquark pairs; (e) diagrams renormalizing the initial-state parton densities. 
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gluons (b), one encounters box diagrams (c), the sphtting of the final-state gluon into gluon 
and light quark-antiquark pairs, as well as diagrams renormalizing the initial-state parton 
densities (e). The evaluation of these amplitudes has been performed in the Feynman gauge 
and the dimensional regularization scheme has been adopted to calculate the singular parts 
of the amplitudes. The masses of light quarks in Fig.|l|(d,ei) have been neglected while 
the mass parameter of the charm quark has been defined on-shell. I have carried out the 
renormalization program in the extended MS scheme in which the massive particles are 
decoupled smoothly for momenta smaller than the quark mass. The exchange of Coulombic 
gluon quanta in the diagram (Ic) leads to a Coulomb singularity ~ 7r^/2/9ij which can be 
isolated by introducing a small relative quark velocity {Sr. Following the standard path 
||18|| , this effect has to be interpreted as the Sommerfeld rescattering correction which can 
effectively be mapped into the cc wave function. As expected, the infrared singularities cancel 
when the emission of soft and coUinear final-state gluons and light quarks, characterized 
by a cut-off A |jl5|,|l9|,^ , is added to the virtual corrections. The collinear initial-state 



singularities can be absorbed, as usual, into the renormalization of the parton densities pT 
defined in the MS factorization scheme. 

The perturbative expansion of the photon-parton cross section can be expressed in terms 
of scaling functions. 
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i = QiliQ denoting the parton targets. For the sake of simplicity, I have identified the 
renormalization scale with the factorization scale ii\ = = . The scaling functions 
depend on the energy variable rj = s/Am^ — 1. c^^ is the lowest-order contribution which 
scales ~ r]~^ ~ Ami/ s asymptotically, c'^^ can be decomposed into a "virtual + soft" (V+S) 
piece and a "hard" (H) gluon-radiation piece. The In-' A singularities of the (V-f-S) cross 
section are mapped into (H), cancelling the equivalent logarithms in this contribution so 
that the limit A ^ can safely be carried out. The nomenclature "hard" and "virtual + 
soft" is therefore a matter of definition, and negative values of c'-^-' may occur in some regions 
of the parameter space. Up to this order, the wave-function at the origin is related to the 
leptonic J / width by 



1 



TT 



with only transverse gluon corrections taken into account explicitly . 

The scaling functions c^i{f]) are shown in Figs.0a/b for the parton cross sections integrated 
over z < Zi where I have chosen zi = 0.9 as discussed before. [Note that the definition of z is 
the same at the nucleon and parton level since the momentum fraction ^ of the partons cancels 
in the ratio z = p ■ k-^ / p ■ k^.] In the range 0.2 ^ rj ^ 2 the hard gluon-radiation piece c^g^^ 
as well as c^^ differ from the curves in Ref. by a few percent since the experimental cut 
z < 0.9 was not implemented properly in one term of 

The following comments can be inferred from the figures, (i) The form of the scaling 
functions resembles the scaling functions in open-charm photoproduction 



However, 

there is an important difference. The "virtual -|- soft" contribution for J/ip production is 
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Figure 2: (a) Coefficients of the QCD corrected total inelastic [z < 0.9] cross section 
•y + g^J/t/j + X in the physically relevant range of the scaling variable r] = s^p/im^ — 1; 
and (b) for 7 + q/q — >• J/ip + X. 
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significantly more negative than for open-cfiarm production. Tlie destructive interference 
with the lowest-order amplitude is not unplausible though, as the momentum transfer of 
virtual gluons has a larger chance [in a quasi-classical approach] to scatter quarks out of 
the small phase-space element centered ed, pc + Pc = VJ/ip than to scatter them from outside 
into this small element, (ii) While c^^^ and c^^^^^^ scale asymptotically ~ 1/s, the hard 
coefficients c^g^-* and c^^* [as well as c^^^] approach plateaus for high energies, built-up by the 
flavor excitation mechanism, (iii) The cross sections on the quark targets are more than one 
order of magnitude smaller than those on the gluon target, (iv) A more detailed presentation 
of the spectra would reveal that the perturbative analysis is not under proper control in the 
limit 2; — i> 1, as anticipated for this singular boundary region |T5[. Outside the diffractive 
region, i.e. in the truly inelastic domain, the perturbation theory is well-behaved however. 



The cross sections for J/ if) photoproduction on nucleons are presented in Figs.^|-p. In 
Fig.^ the leading-order and next-to-leading order calculations are compared with the J/ip 
energy spectra of the two fixed-target photoproduction experiments at photon energies near 

= 100 GeV, corresponding to an invariant energy of about 
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14 GeV. The GRV 
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Figure 3: Energy spectrum da/dz, at the initial photon energy 
with the photoproduction data H,0]. 



100 GeV compared 



parametrizations of the parton densities [Q have been used. They are particularly suited 
to characterize the magnitude of the radiative corrections properly since they allow one to 
compare the results for the Born cross section folded with leading order parton densities, 
with the cross sections consistently evaluated for parton cross sections and parton densities 
in next-to-leading order. As the average momentum fraction of the partons < ^ >~ 0.1 is 
moderate, the curves are not sensitive to the parametrization in the small-x region. Similar 
to the case of open heavy flavor production [^0|,^, the absolute normalization of the cross 
section shows a strong dependence on the value of the charm quark mass. In the static 
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approximation the choice rric = Mj/^/2 is required for a consistent description of the heavy 
bound state formation. However, a smaller mass value might be appropriate for a reasonable 
description of the charm quark creation in the hard scattering process. In order to demon- 
strate this uncertainty the results are shown for two mass values, rric = Mj/^/2 1.55 GeV 
and nic = 1.4 GeV. The i^'-factor, K = ctnlo/ctlo ~ 1-5, consists of two parts, one due to 
the QCD radiative corrections of the leptonic J/ip width ||22| and a second part due to the 
dynamical QCD corrections ||T^, and does not strongly depend on z. The dependence on 
the renormalization/factorization scale Q is reduced considerably in next-to- leading order. 
While the ratio of the cross sections in leading order for Q = rric : (a/2 mc) : Mj/^ is given by 
1.8 : 1.3 : 1, it is much closer to unity, 0.7 : 1.1 : 1, in the next-to-leading order calculation, 
Fig.^. The cross section runs through a maximum |^ near Q \/2mc with broad width. 
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Figure 4: Dependence of the total cross section 7 + P J/ip + X on the renormaliza- 
tion/factorization scale Q at a photon energy of i?^ = 100 GeV. 

the origin of the stable behaviour in Q. In the BLM scheme Q moves from values below 
rric at low energies up to ~ \plm(. at the HERA energy of i/s^p ~ 100 GeV. In particular 
the value at high energies is significantly larger than the corresponding BLM value for J/-?/' 
decays. The typical kinematical energy scale is not set any more by the small gluon energy 
in the J/ tl) decay but rather by the typical initial-state parton energies. I have adopted the 
scale Q = mc in Fig.^ and subsequently. 

In a systematic expansion one may finally add the relativistic corrections as estimated in 
r^ . Two conclusions can be drawn from the final results presented in Fig.|^. (i) The J/ip 
energy dependence dcr/dz^-y+Af J/ip+X) is adequately accounted for by the color-singlet 
model so that the shape of the gluon distribution in the nucleon can be extracted from J/ ip 
photoproduction data with confidence, (ii) The absolute normalization of the cross section is 
somewhat less certain; this is apparent from the comparison with the photoproduction data. 
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[The situation is worse for electroproduction data |TO|JTT|]. However, allowing for higher- 
twist uncertainties of order {A/rric)^ <, 20% for A; > 1, I conclude that the normalization too 
appears to be under semi-quantitative control. 

In Fig.^ I present the prediction of the cross section for the HERA energy range, again 
for two values of the charm quark mass, rric = Mj/^/2 and = 1.4 GeV, respectively. In 
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Figure 5: Total cross section for inelastic J/ip photoproduction ■y + P — > J/ip + X as a. 
function of the photon-proton center of mass energy in the HERA energy range. 

this high energy range the i^-factor is smaller than at low energies, K = ctnlo/clo ~ 0.75, 
a consequence of the negative dip in the c^^^ scaling function of Fig.^]. Note that the LO 
cross section in Fig.^ has been evaluated by using leading-order expressions for the parton 
distributions [^. When adopting the same set of parton distributions for both LO and NLO 
cross sections, the i^'-factor is close to one, depending in detail on the photon-proton center- 
of-mass energy and the choice of the parton distributions. The results in Fig.^ are shown 
for two values of asiV^rric) = 0.25 and 0.31 which correspond to the la lower and upper 
boundary of the error band in Ref . 1 2% , respectively. Since the cross section depends strongly 
on the QCD coupling, I adopt this measured value, thus allowing for a slight inconsistency to 
the extent that the GRV fits are based on a marginally lower value of a^- For = 
one finds, for z < 0.9, a value of about (7(7 + p — J/ip + X) ^ 18 nb at an invariant 7p 
energy of y/s^p ~ 100 GeV; this value rises to about 30 nb if one chooses rric = 1.4 GeV 
and the larger value 0.31 for the QCD coupling. Inclusion of the relativistic corrections as 
estimated in Ref. [|T^ increases the cross section in the HERA energy range by approximately 
10 %. 

Since the momentum fraction of the partons at HERA energies is small, the cross section 
presented in Fig.|^ is sensitive to the parametrization of the gluon distribution in the small-x 
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Figure 6: The total cross section as a function of the photon-proton center of mass energy 
for different parametrizations of the gluon distribution of the proton. 



region < ^ >~ 0.003. This is demonstrated in Fig.|^. The GRV parametrization adopted 
in Fig.^ leads to an almost linear rise of the cross section with the 7p cm. energy. This 
increase is even more pronounced when using the MRS(G) set and considerably less marked 
for MRS(A') p8|. The MRSDO' set gives rise to a much smaller cross section which 



does not strongly depend on the energy in a wide range 30 GeV < ^/s^ < 200 GeV. 
Since the absolute normalization of the cross section is rather sensitive to the value of a.^ and 
the charm quark mass, the discrimination between different parametrizations of the gluon 
density in the proton has to rely on the shape of the cross section [as a function of the 
photon-proton center of mass energy] rather than the absolute size of the prediction. 

The cross sections for inelastic photoproduction of ip' particles can be obtained from the 
results presented here by replacing the leptonic decay width and multiplying with a phase 
space correction factor, a{-fP ip' X) ^ Vt'jVH^ {Mj,^/M^,f x a(7P ^ J/^ X) ^ 
1/4 X aipfP J/ip X). The production of T bottomonium bound states is suppressed, 
compared with J/ip states, by a factor of about 300 at HERA, a consequence of the smaller 
bottom electric charge and the phase space reduction by the large h mass. 

Conclusion: I have shown in this next-to-leading order perturbative QCD analysis that the 
energy shape of the cross section for J/ip photoproduction is adequately described by the 
color-singlet model. A semi-quantitative understanding has been achieved for the absolute 
normalization of the cross section. Higher-twist effects must be included to further improve 
the quality of the theoretical analysis. The predictions for the HERA energy range provide 
a crucial test for the underlying picture as developed so far in the perturbative QCD sector. 
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